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1 . INTRODUCTION 

The spline-under-tension, developed by Schweikert 14 and Cline 2 , 
introduces a parameter which gives some control on the shape of the spline 
curve. The tension spline involves the use of hyperbolic functions and 
can be considered within a general setting proposed by Pruess 13 , wdiere 
other alternatives are discussed. In particular a rational spline due to 
Spath 15 is considered. 

For parametric representations, Nielson 11 describes a polynomial 
alternative to the spline-under-tension. Here use is made of the additional 
freedom given by relaxing C 2 parametric continuity to that of ’geometric’ 
or ’visual’ C 2 continuity. This idea has more recently been taken up by 
Barsky and Beatty 1 . 

The specific problem of shape preserving interpolation has been con- 
sidered by a number of authors. Fritsch and Carlson 6 .and Fritsch and 
Butlana 5 discuss the interpolation of monotonic data using C 1 piecewise 
cubic polynomials. McAllister, Passow and Roulier 9 and Passow and Roulier 12 
consider the problem of interpolating monotonic and convex data. They 
make use of piecewise polynomial Bernstein-Bezier representations and 
introduce additional knots into their schemes. In particular McAllister 
and Roulier 10 describe an algorithm for quadratic spline interpolation. 

In this paper we will discuss a rational spline solution to the problem 
of shape preserving interpolation based on references 3, 4, 7 and 8. The 
rational spline is represented in terms of first derivative values at the 
knots and provides an alternative to the spline-under-tension. The idea 
of making the shape control parameters dependent on the first derivative 
unknowns is then explored. We are then able to preserve the monotonic or 
convex shape of the interpolation data automatically through the solution 
of the resulting non-linear consistency equations of the spline. 



2. A RATIONAL ALTERNATIVE TO THE SPLINE-UNDER-TENSION 


Let (x^,f£),i = 1,...,n, be a given set of real data, where 

x x < x 2 < . . . < and let d^, i = 1,...,n, denote first derivative values 

defined at the knots x., i = 1,...,n. A function s € C 1 [x ,x ] such that 

1 9 n 

(2.1) s(x.) = f. and s^^Cx.) = d. , i = 1 , . . . ,n 

ii i i 99 

is piecewise defined for x £ [x^,x^ +i ], i = 1 , . . . ,n-1 by s(x) = s^(x;r^) , 
where 

(i-e) 2 (i-e+r.e)f .+(i-e) 2 eh.d.-e 2 (i-e)h.d. +e 2 (e+r . ( i -0) ) f . 

(2.2) s . (x ;r . ) = i ± — i ill 

1 1 i + (r.-3)e(i-e) 

and 

(2.3) 0 = (x-x.)/h. , h . = x . -x. , r. > -1 . 

li’i l+l i 9 i 

The parameters r^ will be used to control the shape of the curve s. The 
case r. = 3 is that of piecewise cubic Hermite interpolation. 

The C 2 spline constraints 

(2.4) s ^ 2 ^ (x • ) = s^(x. ) , i = 2 , . . . ,n- 1 , 

l— 1+ 

give the consistency equations 


(2.5) h.d. + [h . (r . -1) +h. (r.-1)]d. + h. d. 

i 1-1 L l 1-1 1-1 l J l 1-1 1+1 

= h.r. A. + h . r. A. , i=2,...,n-1, 

l 1-1 i-i 1-1 l l 


where 


(2.6) A i = (f i+i -f i )/h i , i = 1 , . . . »n-1 

and we assume that d x and d^ are given as end conditions. The case 
Tj. = 3 corresponds to that of cubic spline interpolation. 
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Assume that 


(2.7) >_ r > 2 , i = 1,...,n-1 . 

Then the linear system (2.5) is strictly diagonally dominant and hence 
has a unique solution. The solution is also bounded with respect to 
the r. since 

l 


( 2 . 8 ) 


max |cL| <. [r/(r- 2) ]max | | . 


The rational cubic (2.2) can be written as 

(2.9) s i (x ; r ^) = £ i (x) + e^(x;r.) 
where is the linear interpolant 

(2.10) £.(x) = (1-0)f. + e f. +1 


and 


( 2 . 11 ) 


e d x * r d - 


e(i-e)[(2e-i)(f i+l -f i ) + (i-e)h.d. - eh.d.^] 

1 + (rV-3) 6(1-0) 


Thus, as the parameters r. are increased, it can be shown that the rational 
spline s converges uniformly to a piecewise defined linear interpolant. 

An identical argument applies to the rational spline representation of 
parametric curves . 


Suppose f_^ = f(x^) and f^ = f(x^ +i ), where f £ C 4 [x^x^ +i ]. Then 

an error bound for the rational cubic on [x.,x. 1 is given by 

11+1 ° J 




where 


c. « min{1,(1+r.)/4} and || f || = £ f x I f Cx) | . 

L i 9 i+i J 

This result will influence the choice of the parameters r^ when the 
interpolation data are monotonic or convex. In particular we wish to 
choose r^ such that r^ - 3 = 0(h?) , whilst maintaining monotonicity 
or convexity of the interpolant. 

3. THE INTERPOLATION OF MONOTONIC DATA 
Suppose that 

(3.1) f 1 < f 2 < ... < f (monotonic increasing data) 

and assume the derivative values d^ satisfy the necessary monotonicity 
conditions 


(3.2) 


d^ > 0 , i = 2 , . . . ,n-1 , 


where d, > 0 and d >0 are given. Then 
1 n 


(3.3) 


r . = 1 + (d . + d . )/A. 

l l i+i i 


ensures that s^^xjr.) > 0 on [x.,x. ] and hence the rational cubic 

l l i* i+i 

is monotonic increasing. 

Substituting (3.3) into the C 2 consistency equations (2.5) gives 
the non-linear system 


(3.4) d . [ -c . + a . d. + (a . + a . ) d . + a . d . ]=b. , i = 2 , . . . ,n-1 , 

l l 1-1 1-1 i-i li l i+i l 9 99 9 


where 


(3.5) a. = 1/(h.A.) , b. = A. /h . + A./h. , c. = 1/h. + 1/h. 

l ii* l i-i i-i l l * i i-i l 
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This system has a unique solution satisfying the monotonicity conditions 
(3.2). The "Gauss-Seidel" iteration 


(3.6) 

df k+l) 

1 

= G i (d i-i l) » d i+ ) i ) ’ i=2,...,n-1. 

where 



(3.7) 

G i (C,n) 

’ 2(a. +a.) (t + [t2+4<a i- 1 * a i> b i l!} ' 

1-1 1 

(3.8) 


£ = c . - a . 5 - a .n > 

^ 1 1-1 1 * 


provides a robust algorithm for solving (3.4), being globally convergent 
to the required positive solution. The method also provides an 0(h 4 ) 
accurate solution with reference to the error bound (2.12). 

4. THE INTERPOLATION OF CONVEX DATA 


Let 


(4.1) A x < A 2 < . . . < A ^ (convex data) , 

and assume the derivative values satisfy the necessary convexity conditions 

(4.2) A. < d. < A. , i = 2,...,n-1 , 

l-i ii 

where d. < A., and d > A are given. Then 
i i n n-i a 

(4.3) r £ = 1 + (d i+1 -A.)/(A.-d i ) + (A. - d.) /(d. +1 - A.) 

( 2 ) 

ensures that s. (x;r.) > 0 on [x.,x. ] and hence the rational cubic 

11 1 5 1+1 

is convex. 

The C 2 consistency equations are now 
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giving a non-linear system with a unique solution satisfying the convexity 
conditions (4.2). This solution can be found using a "Gauss-Seidel" 
iteration as in (3.6), where we now take 


(4.5) 


G.(C,n) 


hi (A. -£)*A. 

1-1 1-1 1 

~~ hi (Z - o 

i-i i-i ^ 


+ hi(r ) -A.) i 
1 +hhn -A.) 


A. 

1-1 


I 


and choose intial values 
the residuals alternate, 
monotonic spline is 0(h 4 ) 



Finally, 

accurate 


2,...,n-1, for the iteration such that 
the convex spline method like the 
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